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Abstract - Considering diffusion in the presence of asymmetric disorder, an exact relationship 
between the strength of weak disorder and the electric resistance of the corresponding resistor 
network is revealed, which is valid in arbitrary networks. This implies that the dynamics are 
stable against weak asymmetric disorder if the resistance exponent C, of the network is negative. 
In the case of (^ > 0, numerical analyses of the mean first-passage time r on various fractal lattices 
show that the logarithmic scaling of r with the distance /, Inr ~ Z''', is a general rule, characterized 
by a new dynamical exponent tp of the underlying lattice. 



Introduction. — Disorder — an inevitable feature of 
nature — is known to induce striking slowing down phe- 
nomena in transport processes and in the relaxation of 
systems with many degrees of freedom [T]-[3]. The source 
of complexity in such systems is twofold. First, the dy- 
namics can be regarded as a random walk in the config- 
uration space, which is, in general, rather complicated. 
This problem also arises directly in the context of trans- 
port processes, which usually take place on inhomogeneous 
structures in reality [2]- Second, the system is frequently 
subject to an external source of disorder, such as a ran- 
dom force-field, which can be modeled by quenched (i.e. 
time-independent) random transition rates. The theoret- 
ical description of the dynamics of such systems is a chal- 
lenging task since, in most cases, the complexity induced 
by disorder makes the application of analytical techniques 
of clean systems extremely hard. Diffusion in the presence 
of the most general form of disorder, where transitions be- 
tween pairs of states are non-symmetric, is non-trivial even 
on regular d-dimensional lattices. In the simplest case of 
d = 1, where a single- valued potential exists, the position 
x{t) of the walker varies ultra-slowly with time, obeying 
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(1) 



when the average force acting on the walker is zero [4]. 
Here, (•) denotes an average over different stochastic histo- 
ries, whereas the overbar denotes an average over the ran- 
dom transition rates. This is Sinai's diffusion law, where 



the value tJj ~ 1/2 of the barrier exponent is related to 
the Gaussian fluctuations of the potential landscape. In 
higher dimensions, no potential can be defined, and the 
system is a genuinely non-equilibrium process. For in- 
creasing d, the effects of disorder are expected to be less 
pronounced as there are more and more paths connect- 
ing the initial and final state. Weak disorder expansions 
[5] and perturbative renormalization group analyses [HIT] 
agree in that, for d > dc = 2, the diffusion remains nor- 
mal, i.e. (x2(i)) r^ t, with logarithmic corrections at the 
critical dimension dc ~ 2, while for d < 2 the disorder 
is relevant and results in sub-diffusion. Besides studying 
directly x{t), the inverse question of how large the mean 
first-passage time (MFPT) from a given initial state to a 
final one is frequently asked in such systems [5] . In prob- 
ing the dynamics, the finite-size scaling of the MFPT is 
an alternative possibility as, in general, it obeys the same 
dynamical relation between time and length scales as x{t) 
does. 

In transition networks other than regular lattices neither 
a general criterion for the relevance of asymmetric disor- 
der nor the dynamics in case of relevance are known. Our 
aim in this letter is to study these questions in general. 
We will prove an exact relationship between the strength 
of weak asymmetry in the transition rates and the effec- 
tive resistance of the corresponding resistor network, by 
which a relevance criterion can be formulated in terms of 
the sign of the resistance exponent C. This will be then 
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demonstrated by numerical calculations in various random 
and non-random fractal lattices. In the case of relevance 
(C > 0), the dynamics are found to be logarithmic given 
by Eq.([T]) with a barrier exponent V' that is characteristic 
of the underlying lattice. 

The model and its renormalization. To formu- 
late the above statements precisely, let us assume that 
the system has a finite number of states labeled by the 
integers i = 1,2,...,N and consider a continuous-time 
random walk on them, where the transition rates {pij,pji) 
are from link to link independent random variables. The 
rates Pij and Pji are allowed to be different but their dis- 
tributions are required to be identical, so that the average 
local 'force' is zero, i.e. 



ln{p^j/pj,) 







(2) 



on each link (ij). For the sake of simplicity, we assume, 
furthermore, that every state is reachable from every other 
one. The central quantity in a fixed realization of the tran- 

(A) 

sition network is the MFPT, r^ , which is the expected 
value of the time needed to first reach any state in a fixed 
set A of target sites when starting from state i. These 
quantities with different starting state i obey the follow- 
ing backward master equation [5]: 
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with Ki = 1 for all i and the boundary conditions tI"^' = 
for all i ^ A. As it has been recently pointed out in ref. 
[9] , this form of the equation enables one to calculate t^ 
by recursively eliminating all states other than 1 and A 
one after the other in a way closely related to strong dis- 
order renormalization group (SDRG) methods [3]. When 
eliminating state fc, direct transitions between states that 
were linked with state k are generated with the rates: 



PikPkj/^Pki- (generation) 



(4) 



If the link (ij) has already existed before the elimination 
of k than its rates pij are added to the newly generated 
ones, so that the final renormalized rates read as: 



P,j =Pi3+pI 



(addition) 



(5) 



The parameters Ki at those states that were linked with 
state k are also renormalized as Ki = Ki +PikK}~/ ^iPki- 
When all states except of 1 and A have been eliminated, 
the MFPT can be calculated by 4"^^ = Ki/Y^ieApi^ B 
In the simple case d = 1, the generation rule in Eq. ^ 
is the only operation and its product form results in a 
rapid decrease of effective rates under the elimination (or 



^The denominator here has a direct probabilistic interpretation: 
2j. .pii/2j Pij is the probability that the system, leaving state 
1, will reach state A earlier than state 1. 



renormalization) procedure that signals slow dynamics. In 
more complex networks, also the addition rule given in 
Eq. ([5]) is applied since alternative paths from state i 
to j which do not go through k may exist. The more 
connected the network is the more frequently the addition 
occurs and the slower the effective rate decreases under 
renormalization. 

Renormalization of the weakly asymmetric 
model. — In general, the recursions described above 
cannot be solved analytically. Nevertheless, when merely 
seeking an answer to whether disorder is relevant, the anal- 
ysis can be greatly simplified as follows. It is known in 
the case d = 1 that any weak disorder drives the sys- 
tem to the infinite-randomness fixed point of the SDRG 
[3] that describes the logarithmic behavior in Eq. ([1]). 
Moreover, the behavior of the asymmetry in the effective 
transition rates alone refiects the logarithmic dynamics as 
\n(pii\j/pNi)\ "^ N^/^, and indicates the relevance of weak 
disorder. Therefore, we shall deal only with a weak, asym- 
metric perturbation of the symmetric model of the form 
Pij/ Pji = 1 + eij with eij infinitesimally small, and keep 
track of the renormalization of pij and e^ . The transfor- 
mation rules then read in leading order in eij as 



4 



f-ik + ^kj , 



Pij eij — Pij eij + Pij eij , 



(generation) 
(addition) 



(6) 
(7) 



whereas the rates pij still obey Eqs. (|3][5]), however, they 
are now symmetric, ptj = pji and transform identically to 
the reduction rules of a resistor network with resistances 
rij = 1/pij on links (ij). Let us consider an ensemble of 
networks with fixed resistances r^ and random variables 
eij on each link, for which Eq. ^ implies eij = 0, the 
overbar now denoting the average over e^ on link (ij). 
In order to prove general statements we need to allow eij 
to be non-identically distributed on different links and to 
require that e? = arij , where a is an infinitesimally small 
global constant. 

Relationship with the two-point resistance. 

First, let us consider a special class of networks which 
can be reduced to two fixed states, a and 5, by exclusively 
eliminating states with two links. An example is the hier- 
archical diamond lattice illustrated in Fig. [TJ It is easy to 
see that, here, the generation and addition steps are equiv- 
alent to the well-known reduction of resistors in series and 
in parallel, respectively, and follow the simple rules: 



r = ri+r2, e = ei + £2 



f = r 



(generation) (8) 
^ + 7*2^ , f^ e = Tj^ 61 + r^ £2 (addition). (9) 



When any of the above operations is performed, ei and 
£2 are always independent, £i£2 = 0, so one easily obtains 

that 

efj = ahj (10) 

remains valid at any stage of the renormalization proce- 
dure, all the way to the last link connecting states a and 
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6: 

elb = afab- (11) 

For arbitrary networks, the ehmination of states with 
more than two hnks cannot be avoided in general. When 
decimating a state with n > 2 hnks, the relation Eq. pUj) 
will be broken for the modified n(n— l)/2 links, moreover, 
their rates will become correlated with each other. Sur- 
prisingly, when the network is reduced to two (arbitrary) 
states, a and 5, the relationship in Eq. (fTT|) will still hold. 
In the following, it will be convenient to introduce the 
scaled asymmetry parameters Wy- = Pyfy, which satisfy 
initially wf- — apij- It is sufficient to prove Eq. ([TT|) for 
complete networks (in which every state is linked with ev- 
ery other one) and for an arbitrary set of rates {py }, then 
it applies to any other network by setting pij — a;| = 
for appropriate links, which amounts to deleting that link. 
For A^ = 3, the statement is true since this network be- 
longs to the special class, while for A^ > 3 we shall prove 
it by induction. Assume the Eq. PT|) holds for a fixed 
TV. It follows from Eqs. (jBHT]) that the final asymme- 
try parameter is a linear combination of the initial ones: 

'^ifa (tejl) = '^Kj ^ij i{Pi3})'^ij y where the coefficients 
are functions of the set of rates and the summation goes 
over all links. Using this, Eq. pT|) can be rewritten as 

E^<J[C^^Pi{P^J})?P^, = P^I\{Pv})- Now let US CXtCud 

the network to a one state larger one and keep the nota- 
tion of parameters on links of the A^-state subgraph, while 
denoting the rates and the scaled asymmetry parameters 
on the links from the old state i to the new one by qi and 
(pi, i = 1,2, ... ,N, respectively. If the new state is elimi- 
nated, then effective rates _Py = Qi'lj/J2i=i 1i ^"^^ scaled 
asymmetry parameters Vtij — {qj4>i — <li4'j)/^i=i1i s-re 
generated according to Eqs. dU and ([5]) at all old links, 
and are added to the old parameters as given in Eqs. (O 
and (I7|) . This yields for the parameters of the system with 
N + 1 states: 



^T'\{p^o}. {*}) = E.<, c.M^, + a,,), (13) 

where we have used the shorthand notation Cij for 

Wi \{Pij + Pij})- It is to be shown that the expected 
value of the square of Eq. P^ is equal to Eq. P^ (mul- 
tiplied by a) for any {Pij},{qi\- Using that {wy} and 
{J7ij} are independent and wf- = apij, we obtain that this 
is equivalent to 



«E,<,c^^. = (E.<,c.,%)^ 



(14) 



1 = aqi, the 



Expanding flij in terms of 0^ and using 

r.h.s. assumes the form aX)i=i G^Qi with the coefficients 

Gi = Y^j^i 1jCij/Y,j 9j> where C„ for i > j is defined as 



Cji. Then Eq. ([14)) can be rewritten after some 



algebra in the form 



y^ qiqjqkiCt-j + Cjk + Ckt)"^ == 0, 

i<j<k 



(15) 



where the summation goes over all triangles of links. A 
sufficient condition for Eq. ([T5|) to hold is that 



cl^Hfe,}) + Cf,\{p.,,]) + cif )(fe,}) 



■i] w-t'i] i I '^" ^ jk Wfi^H^^ki Wfiji)—^ (1") 

for all triangles i < j < k (and, in fact, for all closed paths) 
and for any set {pij}- Up to iV = 4 this can be justified 
by direct calculations. The key observation for the induc- 
tion is that when extending the network by adding a new 
state then the coefficients of the old links in the enlarged 
network are modified only by a shift in their argument as 
Clp{{p^J}) ^ CiP{{p,,+P,,}), see Eq. (O. Therefore 
the identity will hold in the enlarged network for all closed 
paths consisting of old links. The same network with iV-|-l 
states can, however, be built from another starting set of 
N states by adding one extra state, as well, in which case 
the identity p^ can be extended to another set of links. 
It is easy to see that applying this reasoning to all possi- 
ble sets of N states (note that a and b are always part of 
this set), the old links cover all possible triangles of the 
extended network and the identity will be valid for iV + 1. 
So, we have proved that Eq. (jlip holds for arbitrary net- 
works. This means that, although the two-point resistance 
depends on the particular network, the variance of the ef- 
fective asymmetry parameter, which can be interpreted as 
the strength of disorder, depends on the resistance uni- 
versally at least for weak disorder. If the resistance is in- 
creasing(decreasing) with the system size then the effective 
asymmetry is getting stronger (weaker) on larger scales. 
So, the dynamical relation between time and length scale 
on lattices with a decreasing resistance is expected to be 
stable against weak disorder; otherwise, it may be altered 
compared to that of the homogeneous system. Measuring 
the distance / between a and b in terms of the linear size 
if the network can be embedded in a d-dimensional lat- 
tice and in terms of the length of shortest path otherwise, 
the resistance scales in many cases as rab{l) ~ I'', where ^ 
is the resistance exponent. The strength of disorder then 
scales with / in leading order as 



\\n{Pab/Pba)\ = ICafcl 



'rab 



(17) 



For C > this is the same type of logarithmic behavior as 
known for d = 1, with the barrier exponent '0'' = C/2. But 
unlike for d = 1, this law will, in general, be modified for 
finite disorder by higher order correction terms which are 
expected to be non-universal. Indeed, numerical results 
that will be presented in the rest of this work indicate, 
in the case of C > 0, a logarithmic scaling law but with 

The results obtained here are in accordance with pre- 
vious ones on regular lattices. For d > 2, the two-point 
resistance tends to a constant as fab{l) ^ l^ + const, with 
( = 2 — d < 0; a.t d = 2, where the diffusion is normal with 
a logarithmic correction, it still increases but just barely 
(formally C = 0), while ioi d = 1, ( = 1. So, in general, we 
expect disorder to be relevant (irrelevant) if C > 0(C < 0) 
while, in the marginal case C = 0, it will induce at most 
corrections to the clean behavior. 
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Numerical analysis. — We have tested this crite- 
rion on fractal lattices, where, in the case of homoge- 
neous rates, the law of diffusion is, in general, of the form 
(x^(i)) ~ i2/d„^ characterized by an anomalous diffusion 
exponent, d^ > 2 ^. This problem has been thoroughly 
investigated but, apart from early Monte Carlo simulation 
studies [TU1[TT| . not in the presence of asymmetric disorder. 
The Einstein relation between diffusion and conduction 
implies that C = d.^ — df, where df is the fractal dimen- 
sion. Prom this, one can see that C > for all fractals in 
d = 2 and, in fact, C is positive for most fractals. 

We have numerically calculated the MFPT by recur- 
sively solving Eq. ([3]) on three different fractal lattices 
with C > 0, namely, the Sierpinski triangle (ST) shown 
in Fig. [T] and bond percolation clusters in 2 and 3 di- 
mensions at the percolation threshold. In the former case, 
the MFPT from one of the tips of the largest triangle 
to the other two tips was considered, whereas in the lat- 
ter case, percolating clusters that connect opposite d — 1- 
dimensional surfaces of a d-dimensional cube have been 
generated (with periodic boundaries in the other direc- 
tions), and the MFPT from a given site on one surface to 
the opposite surface was considered. The rates on each 





Fig. 1: Iteration rule and the 4th generation of the Sierpin- 
ski (a) and the hierarchical diamond lattice (b). The start- 
ing(target) sites are denoted by green(red) dots. 

link were pij = 1, pji = A < 1 with a random (cquiproba- 
ble) orientation. The distributions of the logarithm of the 
MFPT determined in 10^ realizations for different linear 
system sizes L arc found to broaden with L in all three 
cases, see Figs. [212] and to obey the scaling law 



p[Ht/to),L] - L-^p[\n{T/To)L-^] 



(18) 



where the barrier exponents ip are independent of the 
strength of disorder A, and are estimated for the ST and 
d ~ 2 and d = 3 percolation clusters to be, in order, 
0.296(2), 0.46(2) and 0.63(1). 

We have the most accurate estimate for the ST, due to 
an efficient renormalization, essentially by reversing the 
construction procedure. Here, the barrier exponent V' = 
0.296(2) can be seen to be significantly different from -00 = 
C/2==ln(5/3)/ln4R:! 0.368. 

Concerning marginal structures (^ = 0), we have con- 
sidered the hierarchical diamond lattice (see Fig. [1]), 
where df ^ d^ ^ 2 and, in fact. Tab = 1 in each gen- 
eration. The distribution of t/L'^^ is found to converge 
to a limit distribution for L — > oo, but very slowly; the 
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Fig. 2; Scahng plot of the distribution of the logarithm of the 
MFPT on Sierpinski triangles with different linear size L. The 
parameters are tj) — 0.296 and Inro = —1.47. The strength of 
disorder was A = 0.1. Inset: Unsealed distributions. 



typical MFPT, rtyp = exp(lnr), has a logarithmic cor- 
rection of the same form as in d = 2 regular lattices: 
i*"/rtyp ~ D{1 + a/ In rtyp), see Fig. H 

Discussion. — In summary, we have revealed a close 
relationship between the effective strength of asymmetric 
disorder in the diffusion problem and the electric resis- 
tance. This yields a simple relevance criterion in terms of 
the resistance exponent. If C < 0, the dynamics are stable 
against weak asymmetric disorder, while if C > 0, they 
arc unstable and, as numerical results show, the logarith- 
mic scaling is not a peculiarity of one dimension but it is 
the general rule, characterized by the exponent ^p of the 
underlying structure that is independent of other known 
dynamical exponents. In the marginal case, (^ — 0, slow 
corrections to the clean system behavior are expected to 
be general. 

As aforementioned, stochastic processes can be regarded 
as random walks in their configuration spaces. In disor- 
dered stochastic systems with many degrees of freedom, 
logarithmically slow dynamical behavior indeed appears, 
such as in many-particle transport processes with zero- 
range or exclusion interaction |12| . or in the contact pro- 
cess at criticality in one [12] and higher [T3] dimensions. 
The effect of disorder is, however, not restricted to appear 
in one single (critical) point of the parameter space of 
these systems (including the random walk in one dimen- 
sion) but the latter point is surrounded by an extended 
phase, called Griffiths phase [15], where the dynamical 
exponents are finite and vary with the parameters of the 
model non- universally [3]. In the random walk represen- 
tation, the average force acting on the walker becomes 
non-zero in some direction in this phase, therefore this is 
out of the scope of the present treatment (cf. Eq. ([2])). 
Nevertheless, it is plausible to expect systems with a pos- 
itive resistance exponent to have a Griffiths phase, when 
a biased force-field is applied. The non-universal behavior 
with a finite dynamical exponent is, however, suppressed 
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Fig. 3: The same plot for d — 3 percolation clusters as in Fig 
[3 with ii = 0.63 and Inro = -5.7. 



if the lattice contains macroscopic 'dead-ends' (such as 
in critical percolation clusters) from which the walker es- 
capes after an exponentially large trapping time [I], lead- 
ing to logarithmically slow dynamics also in the driven 
phase. The analyses of these phenomena is left for future 
research. 
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